† Corresponding author. E-mail:
Project supported by the Natural Science Foundation of Xinjiang Uygur Autonomous Region, China (Grant Nos. 2018D01C079 and 2018D01C072).
The stable geometries, electronic structures, and magnetic behaviors of the ScLin (n = 2–13) clusters are investigated by using particle swarm optimization (PSO) and density functional theory (DFT). The results show that these clusters have three-dimensional (3D) structures except ScLi2, and ScLi12, and ScLi13 that possess the cage-like structures. In analyses of the average binding energy, second-order difference of energy, and fragmentation energy, ScLi12 cluster is identified as magnetic superatom. The magnetic moment for each of these clusters owns an oscillating curve of different cluster sizes, and their magnetic moments are further investigated using molecular orbitals and jellium model. Of ScLin (n = 2–13) clusters, ScLi12 has the largest spin magnetic moment (3 μB), and molecular orbitals of ScLi12 can be described as
Since the development of cluster science, the physical and chemical properties of small-sized clusters including ScAln (n = 1–8,12),[1]
Recently, attention has been paid to the study of transition metal-doped alkali metals. A strategy to obtain magnetic superatoms was proposed that the magnetic superatoms can be realized by making a transition metal atom locate on the central site of a composite cluster while the surrounding atoms are the alkali atoms.[18–21] A series of transition metal-doped clusters were presented, such as VLin (n = 1–13),[22] ScNan (n = 1–12),[19] VNan (n = 1–12),[18] VCsn (n = 1–12),[18] etc. And the electronic and magnetic properties of these clusters have been investigated by using the first-principles calculations. Previous investigation showed that the electronic states of small metal clusters can be grouped into electronic shell sequence 1S, 1P, 1D, …, much like the way in atoms.[18–20,23,24] Unlike atomic orbitals, the electronic states of superatoms delocalized over the clusters, and the filling of orbitals do not usually exhibit Hund’s rule.[19,20] Reveles et al.[18] found out that magnetic superatoms can be designed by appropriate combinations of localized and delocalized electrons in the valence space of a cluster, and took VNa8, VCs8, and MnAu24(SH)18 as illustrative examples. Zhang et al. pointed out that VLi8 clusters are identified as a magnetic superatom with a jellium sphere electronic structure.[22] Pradhan et al.[25,26] studied Sc-, Ti-, and V-doped Nan clusters in the gradient corrected DFT framework. Pradhan et al.[26] pointed out that ScNa12, ScK12, and ScCs12 clusters each can attain a large spin moment and the addition of alkali atoms to a Sc atom can result in a large spin moment.
It is interesting to note that the first 3d transition-metal atom Sc owns a similar atomic radius to lithium, while to the best of our knowledge there is still no report about the electronic and magnetic structures of Sc–Li binary clusters. Lithium is the first alkali metal in the periodic table whose subshell of valence electrons is 1s22s1. It has been attracting much attention due to its high power and capacity of rechargeable lithium batteries, and the very interesting physical properties of lithium clusters. Zhang et al. pointed out that VLi8 clusters are identified as a magnetic superatom with a jellium sphere electronic structure.[22] Additionally, Li et al. reported that an alternative behavior of odd and even is found in CuLin’s magnetic moment (total or local), for which the clusters with odd number valence electrons show relatively large magnetic effects.[27]
In this paper, the authors investigate in detail the stable geometries and magnetic properties of small-sized ScLin (n = 1–12) clusters. The low energy structures are evaluated by using the particle swarm optimization[28,29] along with the first-principles method. Based on the obtained lowest-energy structures, different growth behaviors and magnetic properties are found. The results show that the core-shell ScLi12 cluster has a regular icosahedral configuration and owns relatively large stability which could be a good candidate for superatoms. The magnetic properties of these clusters can be explained by the jellium sphere model, whose molecular orbitals are in 1S1P1D sequence. The net spin magnetic moment of ScLi12 cluster is 3 μB and its molecular orbital arrangement is
During the calculation, geometric optimization was performed by using the spin-polarized density functional theory (DFT) implemented in the Vienna ab initio simulation package (VASP).[30,31] And the projector augmented wave (PAW) pseudopotentials[32] were used to describe the interaction of valence electrons with the ionic nucleus. The exchange and correlation energies were illustrated by the Perde–Burke–Ernzerhof (PBE) formula of the spin-polarized generalized gradient approximation (GGA).[33] The plane-wave cutoff was established to be 300 eV. The ScLin clusters were placed in a cubic super lattice of 12 × 12 × 12 angstrom (for n < 9) and 15 × 15 × 15 angstrom (for n ≥ 9).
In order to find out the lowest energy structure of the ScLin cluster, we used the CALYPSO[29,34–36] structure prediction method, which is a global minimum search of the free energy surface. It successfully predicted structures for various systems.[37–42] We ran 20 generations for ScLin (n = 2–9) and 30 generations for ScLin (n = 10–13). Each generation produces 20 and 30 structures for n = 2–9 and n = 10–13 respectively. The first generation was randomly generated. For the next generation, 80% of the structure was inherited from the previous generation, and the rest was randomly generated. Then, several structures with the lowest energy were selected as the initial structures, which were re-optimized with more accurate convergence precision. The convergence of the cluster structure was based on the force component less than 10−4 eV/Å and the total energy variation less than 10−5 eV. Subsequent analysis of the cluster structure was performed in the Gaussian09 software package.[43] Harmonic vibrational frequencies were calculated to confirm that the optimized geometry structure is truly minimum and their vibrational frequencies are no imaginary values. The Mulliken population analysis and Adaptive natural density partitioning (AdNDP)[44] bonding analysis were calculated in the Multiwfn program package.[45,46] The AdNDP orbitals were plotted by using the VMD software.[47]
Using the method described above, hundreds of isomers of ScLin (n = 2–13) clusters with different sizes are obtained. After checking in detail the energy of each isomer, the low-energy structures of ScLin (n = 2–13) clusters are obtained. The obtained low-energy ScLin (n = 2–13) clusters are shown in Fig.
As shown in Fig.
The magic number of stable structure is very important for a cluster. Based on the obtained most stable geometry described above, the stabilities of these clusters are further evaluated by using their corresponding average binding energy (Eb), second-order energy difference (Δ2E), and fragmentation energy (Ef), which are defined as follows:
As shown in Fig.
The second-order energy difference Δ2 E values of ScLin clusters are plotted in Fig.
To further probe the stabilities of the ScLin clusters, we calculate the fragmentation energy (Ef) and the results are shown in Fig.
In a word, comparison with the data shown in Fig.
It is interesting to investigate the magnetic properties of transition-metal atom doped alkali clusters such as VLin (n = 1–13),[22] TiNan (n = 1–13),[21] ScNan (n = 4, 5, 6),[25] VNan (n = 4, 5, 6)[25] clusters, etc. Based on the most stable geometries, the magnetic properties of the ScLin (n = 2–13) clusters are investigated. The obtained magnetic moments of ScLin clusters are shown in Fig.
Generally speaking, the orbitals of superatoms (clusters) own complex dependence on not only the valence electron but also the geometry due to the Jahn–Teller effect. The orbitals of superatoms may be explained by using the so-called jellium sphere modelsin which the electron shells can be represented by 1S, 1P, 1D, 2S, … in a similar way to those by atomic orbitals (1s, 1p, 1d,…), and the capital letters are used to separate them from atomic orbitals. In this paper, we investigate in detail the molecular orbitals of the most stable clusters by using the Gaussian 09 code. The obtained orbitals are shown in Fig.
Using the jellium model described above, twe also investigate the molecular orbitals of the other ScLin clusters. For ScLin (n = 2–5) clusters which own about 5–8 valence electrons, the valence electrons are mainly 1S1P states:
To better understand the atomic contribution to the total magnetic moment of each of ScLin clusters, the Mulliken population analysis is further performed by using the Multiwfn software. The obtained atomic magnetic moment is shown in Fig.
The chemical bond and atomic charges are further investigated. Hereafter, we just take ScLi12 for example. The chemical bond, atomic charges of other ScLin clusters are similar to those of the ScLi12 cluster. The chemical bonding analysis of the ScLi12 cluster is performed by using the AdNDP method, which can explain the bonding of a cluster in accordance with n-center two-electron (nc−2e) bonds (1 ≤ n ≤ total number of atoms of the system). As shown in Fig.
In this paper, the stable geometries, the electronic structures and magnetic properties of ScLin (n = 2–13) clusters are investigated by using the PSO method (implemented in the CALYPSO code) and DFT method. The results show that all of the ScLin (n = 2–13) clusters prefer to take 3D structures except ScLi2. Among these ScLin clusters, the ScLi12 and ScLi13 have cage configurations, and the ScLi12 cluster presents a perfect icosahedral structure. After checking in detail the average binding energy, the second-order energy difference, and the fragmentation energy, the ScLi12 cluster is thought to be the magnetic superatom. The magnetic moments of ScLin (n = 2–13) clusters each own an oscillating curve with different cluster sizes. The magnetic behaviors of these clusters can be explained by using the so-called jellium sphere model whose subshells of molecular orbital can be described as 1S1P1D … in the way in isolated atoms. And the ScLi12 has the largest spin magnetic moment (3 μB), and the molecular orbital of ScLi12 can be described as
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] | |
[37] | |
[38] | |
[39] | |
[40] | |
[41] | |
[42] | |
[43] | |
[44] | |
[45] | |
[46] | |
[47] |